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The wetting transition of the Blume-Capel model is studied by a finite-size scaling analysis of 
Lx M lattices where competing boundary fields ±Hi act on the first row or last row of the L rows in 
the strip, respectively. We show that using the appropriate anisotropic version of finite size scaling, 
critical wetting in d = 2 is equivalent to a "bulk" critical phenomenon with exponents a — —1, 
13 = 0, and 7 = 3. These concepts are also verified for the Ising model. For the Blume-Capel model 
it is found that the field strength Hic{T) where critical wetting occurs goes to zero when the bulk 
second-order transition is approached, while H\c{T) stays nonzero in the region where in the bulk 
a first-order transition from the ordered phase, with nonzero spontaneous magnetization, to the 
disordered phase occurs. Interfaces between coexisting phases then show interfacial enrichment of 
a layer of the disordered phase which exhibits in the second order case a finite thickness only. A 
tentative discussion of the scaling behavior of the wetting phase diagram near the tricritical point 
also is given. 

PACS: 68.08.Bc, 68.35.Rh, 64.60F-, 05.70.Fh 



I. INTRODUCTION 

Understanding interfacial phenomena at boundaries [H-Ql such as free surfaces of solids, or walls confining fluids, 
etc., is a topic of very high interest in current condensed matter research, driven by important applications, e.g. in 
micro- or nanofluidics, design of smart nanomaterials etc. But it is also of fundamental interest as a problem 

of statistical mechanics. It is only the latter aspect that will be in our focus here, and hence we shall not deal here 
with wetting properties [9l-[l7| of specific materials, but rather consider wetting in the framework of a simple lattice 
model with nearest neighbor interactions and short-range forces due to the walls. 

In this spirit, a large amount of work has been devoted to the study of simple Ising models with boundary fields; 
see e.g. [l8j for a review, and [l9l - l40| for some pertinent original work. Spin reversal symmetry of this generic model 
implies that phase coexistence between the ordered phases with positive and negative spontaneous magnetization 
occurs at zero bulk field (11=0) in the thermodynamic limit in the bulk; the same holds still true when one considers 
a thin film of a finite thickness L with surfaces at which boundary fields H_\, act, provided one considers the 
antisymmetric situation Hi = —Hl, avoiding "capillary condensation" OHO, HI]) i-e- a shift of bulk coexistence to 
nonzero value of the bulk field due to boundary effects. The latter phenomena is interesting in its own right, e.g. plj . 
and can also be studied in Ising models choosing boundary fields Hi — Hl, e.g. (42| - |43 |. but this problem will not be 
followed up here. 

In principle, the study of a wetting transition, where the boundary induces the formation of a macroscopically 
thick film of the phase that it prefers, coexisting with the other phase separated by an interface at a macroscopic 
distance from the boundary, requires to consider the limit L — !■ oo. For any finite L in the antisymmetric thin film 
geometry, or thin strip geometry, when one considers d — 2 dimensions rather than d = 3, one does not find a wetting 
transition, but rather the " interface localization-delocalization transition" [3 H^-Hlj I37l - l39j : for temperatures above 
Tc{L,Hi) but below the transition temperature Tct of the bulk Ising model, the interface between the coexisting 
phases is "delocalized" , freely fluctuating in the center of the film (or strip, respectively). Here, by center we mean 
the plane (or line) (L -I- l)/2, if we label the layers (rows) of the Ising lattice parallel to the boundaries from n = 1 
to n — L. However, for T < Tc{L,Hi) the interface is tightly bound, and hence "localized", near one of the two 
boundaries. One predicts, however, that for large L the transition temperature Tc{L,Hi) converges rapidly to the 
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wetting transition temperature T,JH^ ) of the semi-infinite system [2^, [30] , and this prediction is consistent with the 
numerical simulations [271 [28l . IslHSSl Issl . |39| . 

One aspect of wetting in the Ising model that found great attention is the behavior near bulk criticality. One can 
show that for a second-order wetting transition the inverse function Hic{T) of Tyj{Hi) behaves as [2l| 

H,,{T)^{T,,-T)^' , (1) 

where Ai is the critical exponent that controls the scaling behavior with the surface field Hi near bulk criticality 
[45.-47.] ■ While Ai in d = 3 is only known approximately from numerical work, in d = 2 Abraham's exact solution 
for the Ising model [l^ implies Ai = i. Note that on the square lattice with exchange constant J between nearest 
neighbors the wetting transition occurs for Hic given as the solution of 

exp(2J/fcBT)[cosh(2J/fcBr) - cosh{2Hia/ ksT)] 

= sinh(2 J/fcsT) , (2) 

and the bulk critical point occurs at exp{2J/kBTcb) = V2 + 1 fisj . 

An intriguing question is, of course, to study how the wetting behavior gets modified if the bulk phase transition 
changes its character. E.g., a simple extension of the Ising model is the Blume-Capel model [1^, [lO], where each 
lattice site (i) can be in one of three states Si = +1, Si = 0, Si = —1 and the Hamiltonian becomes 

n^-jY,s,s, + Dj2s! , (3) 

where the "crystal field" D controls the density of "vacancies" , i.e. sites with Si — 0. 

When D — —00 vacancies are excluded, and the model reduces to the standard two state Ising model, but when 
D / J is of order unity, one has a nontrivial bulk phase diagram with a tricritical point [51] at D — Df. For the square 
lattice, this tricritical point occurs at Dt/J=l. 965(5) and fcsTt/J=0.609 (5). For D > Dt the transition at Tc{D) 
is first order: two phases with opposite nonzero magnetization and the disordered phase (with zero magnetization) 
coexist. 

From the universality principle, we expect that Eq. ([T} with Ai = i would also hold in the second-order region 
of the Blume-Capel model for D < Dt, but to our knowledge this assumption has not yet been tested. Providing 
evidence for this assumption is one of the goals of the present paper. Right at the tricritical point, however, the 
behavior should be different: instead of Eq. ([T]) we now expect 

ific(A,T)cx (Tt-T)^" , (4) 

where the exponent An in d = 3 is expected to be Ai^ = 1/4 [s^. the prediction due to Landau mean-field theory, 
which is expected to hold in d = 3, apart from logarithmic corrections [51[. While in d = 2 the tricritical exponents 
in the bulk arc not those of mean-field theory [5l|, for the Ising model universality class to which the tricritial point 
of the Blume-Capel model belongs fET], they are exactly known from conformal invariance [sst - fssf . E.g. in standard 
notation of critical exponents [56, 57]) at = 8/9, f3t = 1/24, 74 = 37/36, i^t = 5/9, etc. To our knowledge, A^ for this 
problem is not yet known, however. 

When one considers the wetting transition for the part of the transition line Tcb{D) that is first order [D > Dt), 
the wetting transition Hic{T) ends at this line at a nonzero value of the surface field, 

HUD) = HUD,T = T,b{D)) . (5) 
We expect that another power law will exist when D approaches Dt from below, 

HUD)^{D~Dtf , (6) 

but again we are not aware of predictions relating to this critical exponent C,. The evaluation of C, from numerical 
work is a very challenging task, because it requires to locate both the tricritical point and the wetting transitions 
with very high accuracy. 
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Finally, a very interesting aspect relates to the character of the interface between the coexisting phases with 
positive and negative magnetization, when one approaches the bulk pha se transition: then interfacial adsorption of 
the third phase (the disordered phase) at the interface can occur [58l - t60| . If we denote the coexisting bulk phases as 
A and B, approaching the bulk transition line in the second-order region one expects an interfacial wetting transition 
A\B — >■ A\DO\B, where DO stands symbolically for an intruding layer of the disordered phase having predominantly 
states where Si — dominates at the interface. To quantify this effect one considers the net adsorption defined by 

M 

Wo{T)^N-^Y.^{6o,s,)i\-i-{So,s.)i\i] , (7) 

i=l 

where N is the total number of lattice sites, M is the number of lattice sites per row parallel to the walls, S^^p is the 
Kronecker symbol, and (• • • )i|_i a statistical average in the presence of an interface {A — 1, B = —1), while (• • • )i|i is 
the corresponding average for an equivalent system but without interfaces. Since 150,5; — ^ — Sf, one can conclude also 
that Wo{T) = (daAs/dD) where ctab is the interfacial tension for the considered interface. Since for D < Dt the bulk 
phase transition of the Blume-Capel model falls in the Ising universality class, we know that aAB oc (1 — T/Tcb{D))'^ 
with fi — {d— = 1 for the Blume-Capel model. Writing t = 1 — T/Tcb{D), we find 

thus if we define a divergence of interfacial adsorption via 

WoiT)^t-^ , (9) 

we find oj = for D < Dt, i.e. in the Ising-like regime there should be no critical divergence of the net adsorption. 
However, right at the tricritical point, the behavior should be different: from v = 5/9 as quoted above we conclude 
Wo{T) cx with UJ = 4/9 ioT D = Dt. This result is consistent with Monte Carlo results of Selke et al. [59l |. 

For Dt = 1.965 < D < 2 the transition of the Blume-Capel model is believed to be first order, and the emerging 
interfaces A\DO and DO\B remain sharp as T ^ Tc{D). Reducing the problem to interfaces in the solid-on-solid 
model, one can argue that Eq. © also holds in the first-order region, but ui — 1/3 fl2j. These considersations have 
been checked by early Monte Carlo work [58r,60j for interfaces between coexisting bulk phases, but no study in the 
context of wetting at external boundaries has as yet been performed, to our knowledge. 

In the present work, we wish to contribute filling this gap and study wetting behavior for the two-dimensional 
Blume-Capel model in the case of thin films for which surface fields at the boundaries act, as described in Sec. 2. 
However, in order to do so, we have found it necessary to reconsider the simulation methodology for the study of 
critical wetting. In fact, by applying Monte Carlo simulations, we shall necessarily study finite systems, namely strips 
of width L (in y-direction) and length M (in x-direction, where a periodic boundary condition acts) . As a consequence, 
finite size effects matter, and hence in the next section we also shall give the background on the proper finite size 
scaling analysis of such simulation "data" , for the case of critical wetting in d = 2 dimensions, taking into account 
that one deals there with an anisotropic critical phenomenon of a special character: if we consider the magnetization 
m of the strip as the "order parameter" of the transition, its critical exponent f3 = 0. This fact does not seem to 
have found much attention in the previous Monte Carlo studies of critical wetting in d = 2 [2^, [s^ . This revised 
methodology for the study of critical wetting, which is outlined in Sec. 2, is another important result of our work. 

Sec. 3 then presents our Monte Carlo results, which then are analyzed according to these concepts, and the wetting 
transitions are discussed for a variety of choices of D, both cases D < Dt, D ^ Dt and D > Dt will be considered. 
Sec. 4 then discusses our result on interfacial adsorption, and Sec. 5 summarizes our conclusions. 

II. THEORETICAL BACKGROUND AND DETAILS ON SIMULATION METHODS 

A. The model 

We consider the Hamiltonian of the 3-state Blume-Capel model [i^, [13] where each lattice site i carries a spin Si 
that can take on the values 5*^ = ±1, 0, at a square lattice in a L x M geometry, where periodic boundary conditions 
act in the x-direction (where the lattice is M rows long), while free boundary conditions are used in the y-direction, 
where boundary fields Hi, act on the first and the last row. Thus the Hamiltonian is 
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i iG rowl rowL 

Here J is the exchange constant between spins at nearest neighbor sites, which we take homogeneous throughout 
the system. Thus, we disregard the possibihty to take the exchange in the boundary rows different from the 
exchange J in the interior of the system, that is frequently considered in studies of wetting phenomena in the Ising 
model [18i, ,25,. 26. 31 33. 39], The bulk field H acts on all lattice sites, while the "surface fields" [45-47] only act on 
the spins in the first (1) and last (L) row, where the free boundary conditions apply. 

We also specialize on the particular antisymmetric situation Hi — —H^ < and consider the thermodynamic 
limit {L oo, AI — oo) before we consider the limit ^ 0+. Then, the system undergoes two phase transitions: 
at the temperature Tcb{D) the phase transition occurs from the disordered "paramagnetic" phase to the ordered 
"ferromagnetic" phase, where we have used the terminology of the Ising model. Taking the lattice spacing as our unit 
of length, the total number of spins is just the product of the linear dimensions of the system, N — LAI. In fact, 
when we define the magnetization m per lattice site as 



^ = ^T.s^ . (11) 

1=1 

its thermal expectation value (to)t for temperatures T < Tcb{D) will be nonzero and positive in the considered limit. 
However, a second phase transition occurs at a lower temperature T„(i7i) for small enough absolute values \Hi\ of 
the surface field: Note that the surface field Hi is oppositely oriented to the positive bulk field H, but — > 0+ while 
Hi, stays finite. Thus, near Tcb{D) the surface field stabilizes a macroscopically thick layer of negative magnetization 
near the lower boundary, where Hi < acts, separated by an interface from the bulk, where the magnetization is 
positive. At r„(iJi), a transition occurs where this interface gets localized near the lower boundary: in the extreme 
case, the domain with negative magnetization disappears completely. In the Ising model, which results from Eq. (|10p 
as the limiting case D -> — oo, this wetting transition is second order throughout the regime < \Hi\ < J. 



B. Critical wetting in d — 2 dimensions: A brief review 

For the semi-infinite system described above, the wetting transition is a singularity of the surface excess free energy 
fs^\T; H, Hi), defined from standard decomposition of the total free energy into the bulk term and boundary terms, 
for L — > cx), 7\f — > oo 

FiT, H, Hi,Hl,L, M)/{LM) = U{T, H) + ^fi^\T, H, Hi) 

+ ifW{T,H,HL) . (12) 

The singular part of this boundary free energy, also called "wall tension" or "wall excess free energy" [l|, 0, 19l-[l3. [l8| , 
near T^{Hi) is expected to satisfy a scaling behavior [H, HH [13, [Ml Is!] where i = 1 — T/T^{Hi) — > 0, 

/ilg/fcsT=Kp-"=F,(i7|r^0 • (13) 

Fs being a scaling function that we do not specify here. This is analogous to standard scaling in the bulk near T^b 
[3, where r = 1 - T/Tcb ^ 0, 



/Ming/fcBr^|r|^-"F,(i/|rr^) . (14) 

In Eq. ([13]), a is the specific heat exponent, and the "gap exponent" A is related to the standard exponents (3 of 
the spontaneous magnetization {{rnjT oc r^) and 7 of the susceptibility {x = {d{m) /dH)T) oc \t\~'* via the scaling 
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relation A = 7 + /3 (Sg. Thus as, are analogous exponents characterizing the singular behavior of the critical 
wetting transition in d = 2 dimensions. 

As always when one considers a critical phenomenon, a diverging correlation length exists. Since wetting can be 
viewed as an "interface unbinding" transition [qL^^ when approached from below, it is natural to study the correlation 
function C{x) describing the correlation of fluctuations of the height of the contour £{x) separating the domain with 
negative magnetization near the lower boundary from the region with positive magnetization in the bulk. Defining 
S£{x) as Si{x) = i{x) - {£), this yields 

G{x) = {Se{Q)6e{x))T , (15) 

noting translational invariance in a;-direction. By (...)t we denote an average in the canonical ensemble at temperature 
T; in the following the subscript T will be omitted. Near critical wetting in d = 2 G{x) takes the scaling form j61|J 

G(a;)-x-(i+''N)g(:,/^ll) , (16) 

writing analogously to the decay of bulk correlations near criticality, G(r) = '■'^^^^'''^(j'/^fc); with 77, 771 1 being 
appropriate exponents describing the variation of these correlations at large distances right at bulk criticality or 
critical wetting, where (or ^||) is infinite. For ^||, a scaling relation analogous to Eq. (IT51) holds, 

en =t-''N^~|(Ht-^.) ; (17) 

again this Ansatz is inspired by the corresponding bulk behavior, 

^b = r-''UH\T\-^) , (18) 

and again are scaling functions that we do not specify here. 

We disregard here the problem that in general a nonlocal theory [63l - [6^ is required (i^l , since it does not change the 
following conclusions. Capillary wave theory implies that 7711 =0, so there is no second independent critical exponent 
as in the bulk. For a more detailed reasoning to explain why for the case of critical wetting with short-range forces 
the exponent 771 1 is zero in all dimensions, and hence there is a single independent critical exponent, we refer to the 
literature [l^]- Remember that using hyperscaling relations dv = 2 — a — ^ + 2/3 and the scaling law 7 — v{2 — 77) 
one can express all static critical exponents just in terms of only two independent exponents v and 77 in the bulk [56|. 
For critical wetting, in d = 2, there is only one independent critical exponent namely, [ll| 

^11=2 (19) 

and the hyperscaling relation for interfacial phenomena in d — 2 dimensions v\\ — 2 ~ oig [l2| then implies — 0. 
From Eq. p6|) we then conclude that on a length scale x the interface exhibits (at t = 0) a mean-square displacement 
also proportional to x, while for t > one concludes that on the length scale the mean-square displacement in 
y-direction also is proportional to ^n, while for x S> the correlation G{x) — > 0. This leads to the conclusion that 
the correlation length describing the fluctuations of the interface in y direction scales as 

C_L cx i-'^^(fori7 = 0),i/_L = 1 ■ (20) 

Noting from Eq. (fT3|) that surface excess magnetization 771^ and surface excess susceptibility Xs follow from deriva- 
tives with respect to the field 

m, = -dfi'HT,H,H,)/dH\T,H, , (21) 
X. = -id-'fi'\T,H,H,)/dH')T,H, (22) 
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we conclude that these quantities have the singularities 



m, cx t 



2-a.-A. 



m, oc 



(23) 



and 



X. cxt^-"=-^^= • Xscxi-^' . (24) 

Using then the consideration that this excess susceptibility is just caused by displacements of the interface between 
the domain of negative magnetization at the boundary and the bulk, one can use Eq. (jl6p to deduce a further scaling 
relation between , i/| | and r]\ | , namely [13] 

A, = (i.||/2)[(d-l) + 2-77||]=3, (d = 2) . (25) 
Also the excess magnetization obviously simply is related to and we have, in d = 2, 



nis oc {£) oc oc ■ (26) 

Thus, near to critical wetting in d = 2 the typical mean distance of the interface from the wall and the typical 
excursions of this interface from its mean are of the same order, iy± = —/3s = 1. 



C. Finite-size scaling: wetting vs bulk transitions 

All above relations did refer to the case that we were considering the limits Af — > oo and L ^ oo first, and then 
consider the limits — >■ 0, < — >■ near the wetting transition. While such an approach is natural in the context of 
analytic theories (lol - [l2l |. it it not sufficient to understand simulations, where we wish to take these limits in reverse 
order. When L is kept large but finite, it is clear that for temperatures above the wetting transition and H = 
there is no physical distinction between the domain with negative magnetization in the lower half of the system 
and the domain with positive magnetization in the upper half: Previously, for the semi-infinite case, the positive 
magnetization in the bulk was singled out by taking the limit i7 — > 0+ first, and then the positive boundary field Hl 
for L — >■ CX) was irrelevant, while for L finite it is critical to maintain the symmetry of the situation with respect to 
the sign of the magnetization ioi H = 0. As a consequence, we conclude that for the present situation there is no 
bulk magnetization to consider, and it rather is the total magnetization of the system or strictly speaking, its absolute 
value, which undergoes a transition from zero to a nonzero value when we cross the wetting transition temperature, 
and extrapolate simulation results towards i — )■ oo. 

As a consequence of this consideration, we propose a scaling assumption for the distribution function Pl.a/ (?7i) of 
the total magnetization in this finite geometry as follows [67, . ,6S] 



PLM^n)=^p^P{L'^ii/-^/M,M/C\\,m^'P^) . (27) 

Note that this expression generalizes the standard expression for finite-size scaling in isotropic systems that have 
linear dimension L in all spatial directions [7l|, as appropriate for phase transitions in the bulk 



PL(™) = cf/^n(V6,™4^/') (28) 

to systems with anisotropic linear dimensions L,M and anisotropic correlation length exponents While for 

isotropic critical phenomena an anisotropic system shape would lead to a dependence simply on the "aspect ratio" 
L/M of the system, the fact that M scales with ^|| and L scales with can be used to demonstrate [13,1111 that the L- 
dependence enters via "the generalized aspect ratio" L'^n/'^^/A/ in the scaling function, Eq. (|27t . With — 2, v± = 1 
we simply have an argument L'^/M, which needs to be kept constant when the variation with M is studied. Such 
gen eralizations of finite-size scaling to anisotropic criticality have been discussed earlier for the Kasteleyn transition 
|69l [70| . which also exhibits = 2u± and hence the same generalized aspect ratio /AI applies. The prefactor 
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^[j'^''" in Eq. ([27]), and likewise the prefactor ^1''' in Eq. ([28]), ensures that the probabihty Ph.uijn) can be properly 
normalized, 

+1 

PL,Mim)dm:^l . (29) 
-1 

Taking now suitable moments of Pl,a/(?7i), analogous to the procedure at isotropic phase transitions in the bulk 
[ll III, we obtain 

+1 

(|to|) = J dm\m\PL^M{m) = 
-1 

= CiT'^^''' ' m(L''i I /''^ /M, M/C\ I ) , (30) 



fc = l,2,... . (31) 

In particular, from Eqs. (|5T|) we derive the standard behavior for the "susceptibility". We denote it by x' 

since rather than (m)^ is subtracted, see Ref. [zlj for a discussion. So, one has 

fcsTx' = iAf((m2) - (|m|)2) , (32) 

and hence 

knTx' = LM^;^'^^''"x{L''u/'^^/M,M/^\\) , (33) 

with the scaling function x = — (?ti)^, omitting all arguments for simplicity. 

Now the key task is to identify the "order parameter exponent" /3 in the context of this description of a critical 
wetting transition. For this purpose, we note that the singular behavior of x' can only be due to the singular behavior 
of Xs oc t~'^, Eq. ([M)) . where we used the results = 0, A = 3. A key point to note is that x' was normalitzed per 
spin, relating to the total volume LM, while Xs in Eq. (j24p is taken relative to boundary sites, so is only normalized 
by M. Hence we conclude that at Tyj, we have, using L^/M = c in Eq. (1551) to eliminate L 

ksTx' oc v^M3/2-2/3/-ii (34) 
irrespective of how the finite constant c is chosen. Writing then 

ksTx'lT^ = kBTxs/L = kBTxsM-^'^c-^'^ (35) 
and using a finite-size scaling relation for Xs 

Xs^t-^x{Ml£,\\)^ilx{M/^\\)^M^ ,T^T^ (36) 

and hence we find 

fesTxlT™ oc M3/2(cx i^) . (37) 
Comparison of Eqs. ([Ml) and (I57|) yields a central result of this section, namely 
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/3 = , (38) 
and Eq. (I37p also can be interpreted in terms of the standard finite-size scaling result 

^bTxIt^ oc M^/'^ii , 7 = 3. (39) 

Second order transitions with an exponent (3 = are rather unusual; for another recent example see Jaubert et al. 
[73} . In the present case, we can also understand it from the fact that (|to|) must tend to nii, in the thermodynamic 
limit in the partially wet phase, where the interface is bound to one of the walls up to the transition point, while in 
the wet phase for the considered limit (|m|) = 0. It is interesting to note that these critical exponents also satisfy the 
usual scaling relation with the gap exponent As introduced in Eq. (|13p . 



7 + /3 = As=3 , (40) 
as well as the hyperscaling relation for anisotropic bulk critical phenomena in d = 2 dimensions, 

i^ll +i^_L = 2/3 + 7 = 3 ■ (41) 

One should not confuse the exponent /3 of the total magnetization, resulting in the limit i — > cx), M — 00, 
L^/M = c, with the exponent (3s of the surface excess magnetization rUg defined in Eq. ([25)1 : while (m) — ?► as i — 0, 
TTis diverges as t — > 0. Actually, we have the standard scaling relations 

(3s = (^ , -fs = 1 + v± , as = a + v± (42) 

with a = —1, as is easily checked by noting that 7 + 2/3 = 2 — a = 3. 
A short manipulation of Eq. shows that it can be rewritten as 

fcBrx' = i^/''-^(i'^||/''VM,i/e±) , (43) 

where x is another scaling function. For a fixed generalized aspect ratio i'^n/'^^/Af (= c) this function has a maximum 
at same value X„iax of the argument X = L/£^±. This implies that the height of the maximum scales as ksTXmax — 
jji in the present case); see also Eq. ((37|) : and its position Tmax scales as 

T/Tmax - 1 (X L-^/''^ {= L-i). (44) 



Incidentally, this extrapolation has already tentatively been used in an early work 27[ without detailed justification. 

The interpretation of the wetting transition in analogy to bulk critical phenomena, considering a particular way of 
taking the thermodynamic limit (cf. Fig.[T]), is reminiscent of studies of filling transitions in double wedge 74, 75] or 
bipyramid [tI, [tJ geometries; however, there the critical exponents are different. A useful consequence of Eq. ((38|) is 
that both (|m|) and {m?'^) right at the critical wetting transition become completely independent of linear dimensions, 
but depend still on the constant c = i^/M, the generalized aspect ratio. Plotting (|m|) vs. temperature for different 
choices of M at constant c, critical wetting should show up via an unique intersection point. The fact that for wetting 
transitions in d = 2 for the choice M/L^ = const, a unique intersection point of (|m|) versus T curves occurs has been 
only noted previously [zl] in a study where long range boundary fields were applied, where = 00 holds [zll, and 
hence the critical behavior is rather anomalous. Fig. [T] shows a schematic sketch where the typical excursion of 
the interface is still smaller than i, by a factor of about 2.5, and similarly C|| is smaller than M. The absolute value 
(|m|) of the magnetization then still is of the same order as the bulk magnetization mb that one encounters inside of 
the domains. The distribution of ni then is distinctly bimodal, and (m) = 0, since the interface is bound to the upper 
wall with the same probability as the lower wall. However, when T approaches T^,, the attractive force between the 
wall and the interface is compensated by the entropic repulsion, the excursions of the interface extend over the whole 
width L, and then show a correlation over a lateral length M (x L^. This interface configuration at T^{Hi) then is 
selfsimilar, irrespective how large L is chosen, as expected for critical phenomena. Thus, the result (3 = must not 
be confused with a first-order wetting transition, where one would have a discontinuous change from a bound to an 
unbound state of the interface, rather than this continuous unbinding. 
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Of course, it is of interest to ask how the behavior changes when one does not take the thermodynamic hmit in this 
particular way where the generalized aspect ration LP' /M is fixed, but rather choosing one of the linear dimensions 
large but finite, and varying only the magnitude of the other linear dimension. Thus, if L is kept fixed and one 
considers the limit Af — ^ c», the system becomes quasi-one- dimensional. While mean- field theory would predict a 
sharp interface localization-delocalization transition in this limit 30], one expects that thermal fluctuations cause a 
rounding of this transition since the system is equivalent to a one- dimensional Ising system H, and this expectation 
has in fact been verified numerically [79j . If instead we keep M finite and increase L, the growth of ^|| would be 
limited by the magnitude of Af , and hence also as well as the distance of the interface from one of the walls in the 
bound state would be limited by a value proportional to a/M, irrespective of L, while in the wet state its equilibrium 
position is at L/2. As a consequence, data taken in such a way would appear as if the wetting transition were weakly 
of first order, which it is not. Thus, we conclude that only the finite-size scaling analysis at constant generalized 
aspect ratio I? jM = const, is the most useful choice. 



D. Comments on the simulation parameters 



In most of the numerical work, we have chosen a particular value of this generalized aspect ratio, namely 



L'^/M = c = 9/8 . (45) 

Of course, the value of the constant c in principle is arbitrary, and the results on the location of T,^{Hi) 
and the critical exponents should not depend on this choice. Our motivation for this choice was simple that it 
yields solutions for L and M which are both integer, as it must be on the lattice, and rather small: (L, M) = 
(6, 32); (12, 128); (18, 288); (24, 512); (30, 800); and (36, 1152). However, the effect of considering other choices of c will 
be discussed below. 

Monte Carlo simulations were then performed using the standard Metropolis algorithm, see e.g. [zl] for a review. 
Typical runs are performed over a length of 10^ Monte Carlo steps per lattice site (MCS), disregarding the first 2 x 10^ 
MCS to allow the system to reach equilibrium. Note that for systems far below bulk criticality exposed to boundary 



fields cluster algorithms do not present any advantage [44 1 



In addition, part of the bulk phase diagram of the Blume-Capel model has been calculated, by applying again 
standard Monte Carlo methods j72| but using finite-size scaling for L x L lattices with periodic boundary conditions 
in both X and y directions, see Fig. [5J For earlier studies see e.g. [sol. IsTj. 



III. MONTE CARLO RESULTS ON THE WETTING BEHAVIOR 

Motivated by the analysis of Sec. II, we have analyzed the first two moments (|m|), (to^) and the cumulant fnlW'. 



UiT) = l-{m*)/[3{m'f] , (46) 

of the magnetization distribution and plot them versus temperature T for a few typical values oi D/ J, as shown in 
Figs. [3]- [5]) . As expected, rather well-defined intersection points in all three quantities (\'rn\), {m?) and U{T) can be 
found at the same estimate kBTyj{Hi) / J within reasonably small errors. As always [Zll, [23> the statistical accuracy 
is better for (|to|) and (m^) than for the cumulant. For standard phase transitions in the bulk [t^I, and (m^) 

would not exhibit "universal", size-independent, intersection points, of course. So, the behavior seen in Figs. [S]- [S] is 
clear evidence for the scaling description of Sec. II, and also the direct observation of configuration snapshots, see 
e.g. Fig. [51 is compatible with the qualitative picture that was developed (Fig. [T]). Note that in the wet phase the 
average position of the interface is in the middle of the L x M strip, in between the rows y — L/2 and y = 1 + L/2, 
for L even as chosen here; but due to capillary waves on the scale M in x— direction the interface makes excursions 
of order a/M in y— direction, which are of the order of L/2, for our choice of geometry (Eq. (|45I) V 

On the other hand, in the case shown for D = — oo, i.e. the Ising limit of the Blume-Capel model (Fig. [3]), the 
location of the wetting transition is exactly known [l^], Eq. and highlighted by a vertical straight line. Evidently, 
the intersection of the curves in Fig. [3] occurs at a temperature compatible with this prediction. 

We do expect that the critical wetting transition of the Blume-Capel model falls in the same universality class as 
for the Ising model. Thus, for the same choice of the generalized aspect ratio c, Eq. P5|) . we expect the same ordinate 
values of the intersection points of (m^) and U{T), respectively. Within our accuracy, the data are compatible 

with this expectation. 
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Fig. [7] analysis the scaling of the peak heights and peak locations of the susceptibility ksTx', ci. Eq. (p7)) and 
Eq. pi)) . respectively. Also, Figs. [5] and [3] present some examples for the "data collapse" on master curves when 
(|m|), (to^) and U are plotted as a function of t^/M as a test of Eq. ([50]) . (|m|) = rh{L'^/M, Mt^S^^^) where we used 

the asymptotic power law for (^||, ^|| = = ^||^^^- One sees that a fair data collapse does in fact occur, although 

deviations due to both statistical errors and systematic effects are present, when M and/or ^|| are not large enough or 
when data too far away from Tw{Hi) are included. Of course. Figs. [7][9] are just examples only, but representative for 
the general pattern of behavior. In any case, we do conclude that using such finite-size scaling analyses as presented 
here one can locate wetting transitions for two dimensional lattice models, such as the Blume-Capel model on the 
square lattice, with reasonable accuracy. 

On view of the well-defined intersection points and the data collapse in all three quantities (|to|), {rn?) and U{T) 
so far studied, figures [3][S] and [SJHl respectively, it is worth to analyze the scaling behavior of the distribution function 
Pl,m{^) of the total magnetization given by equation (|27p . In fact, since one has /3 — the prefactor and the third 
scaling argument are constants. Also, just at ksT-uj^Hi)/ J (i.e. t = 0) the second scaling argument vanishes due 
to the divergence of the correlation length ^[|, so that the distribution function depends on the generalized aspect 

ratio c — L^w/'^^ /M , equation (|^5|) . of the sample only. Figures HU] (a) and (b) show plots of PL.Aiim) versus m as 
obtained by keeping c = 9/8 and for the cases D/J = —oo and D/J = 1.50, respectively. It might be expected that 
the curves of the probability distribution would be independent of D, except for normalization corrections due to the 
different density of vacancies. Since we found that the shape of the probability distribution depends sensitively on T, 
it could be that the qualitative differences between the curves shown in figures [TU] (a) and (b) would be due to the 
uncertainties in the location of ksT^iHi) / J for the case D/J — 1.50. 

On the other hand, by keeping the temperature, the crystal field, and the surface magnetic field just at the wetting 
transition point, we have measured the dependence of PL.Aiifn) on the generalized aspect ratio, for the Ising model, 
as shown in figure [11] (a). In that figure we choose i = 30 while M is varied. We recall the non-trivial shapes of 
the probability distributions that exhibit a single peak around m = for rather elongated samples (i.e. c < 0.75), 
which monotonically crosses over to a bimodal distribution in the limit of more 'cubic' samples (i.e. c > 1.80). 
Also, the generalized aspect ratio selected for our detailed Monte Carlo simulations (c = 9/8), which can roughly 
be located within the crossover regime, exhibits not only the central peak, but it also shows the onset of growth of 
lateral shoulders close to m ~ ±0.60. By using the distributions already shown in figure [TT] (a) we can also measure 
the dependence of the crossing points on the generalized aspect ratio of the sample (figure [TT] (b)). From this plot it 
follows that the suitable range of c in order to obtain a reliable intersection point for both (|m|) and U{T) in Monte 
Carlo simulations is rather narrow, say 0.75 < c < 2, and slightly broader for the case of (m^). In order to be able 
to locate intersections with reasonable accuracy, it is advisable that the cumulant intersection is neither close to zero 
nor close to its maximun value (2/3). It is also reassuring to see that for no shape the distribution at the wetting 
transition resembles the shape expected for a first-order wetting transition, which would be three delta functions at 
zero and positive and negative spontaneous magnetization in the thermodynamic limit, and a distribution with three 
peaks of approximately Gaussian shape in a large but finite system. 

The location of the wetting transitions as a function of two parameters. Hi/ J and D/J, by determining well- 
defined intersection points still requires extensive computations and hence is a challenging task. Figs. [T^] and 1131 
Since Hic{J) at ksTt/ J = 0.609 decreases only rather slowly with increasing D when one approaches Dt, it obviously 
is very difficult to estimate the exponent Ai^ (Eg. ffl numericallv. In mean- field theory, supposed to be exact for 
d > 4 we know that Ai = 1/2 and Au = 1/4 [ZtI. Is^ls^: while in d = 3 dimensions numerical estimates yielded [2^ 
Ai « 0.45, and we know that in d = 2 Ai = 1/2 again. So, we see that the exponent Ai never deviates much from 
its mean-field value, if at all. If we speculate the same observation to be true for Ait, we would expect that the data 
in Fig. [T3]for T — Tt vary as (1 — D/DtY^'^, which is not unreasonable. 

We also see that in the region where the transition in the bulk is first order, the wetting transition lines Hic{D) 
indeed end at nonzero values H^J^D) at the bulk transition line, and these values decrease as one approaches the 
tricritical point, qualitatively compatible with Eq. (|5]). Thus, we can sketch the global wetting behavior as shown 
schematically in Fig. 1141 however, it would be premature to attempt to estimate the exponent <j: our data are clearly 
too limited for this purpose, and we have not tried to extend our study, since we feel the incomplete knowledge on 
the precise location of the bulk tricritical point [13, [Ml would hamper such a study. 

As a final problem we consider the enrichment of vacancies at the interface between oppositely magnetized domains 
in the Blume-Capel model. This problem has been considered earlier by Selke et al. [58l - l60| for the case of free 
unbound interfaces, while here we consider the extension where the interfaces are confined between competing walls 
and may undergo a wetting (interface unbinding) transition. Fig. [15] shows plots of the profiles of the vacancy density 
(a) and the magnetization (b) for a situation of incomplete wetting. One observes a clear enrichment of the vacancies 
in the interfacial region. The more the interface unbinds from the wall, the more the peak of the vacancy concentration 
moves away from the wall (Figs. [15] a, b). Note that for the low temperature shown, the magnetization in the bulk 
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tends to (m) = — 1, i.e. the vacancy concentration inside the bulk is very smalL 

Close to the wetting transition, the interface gets detached from the wall, and then the interface is located in the 
middle of the strip, at z = (L + l)/2 on average. However, since the interface strongly fluctuates around its average 
position a straightforward measurement of the vacancy concentration profile would yield an almost horizontal flat 
curve across the strip. So we have defined a local coarse-graining of the interface, taking segments of length Ax — 8 
in the x-direction, and determining the local center of mass z = £(x) of the interface in each segment, we compute the 
vacancy distribution in each segment separately, relative to the center in each segment, and superpose the distributions 
from the individual segments such that their centers coincide. In this way one obtains vacancy profiles with a clear 
peak in the center of the strip in the wet phase (Fig. [TSh). Finally, when one reaches the phase boundary of the bulk, 
the disordered phase takes over in the film, and the vacancy concentration becomes (almost) unity in the system, 
except near the boundaries where the surface fields still stabilize layers of up-spins and down-spins, respectively 
(Fig. [TSH). Note, of course, that the details of the curves in Fig. [TSb.d do depend on this coarse-graining length Aa; 
distinctly. 

This ambiguity that results depend on the coarse-graining is avoided when one simply computes the average fraction 
of vacancies in the system (Fig. [T6| . While for |iJi| chosen such that the system is in the incompletely wet state up 
to the transition in the bulk one sees also a clear jump in the vacancy concentration from a small value to almost 
unity when this transition in the bulk occurs, a much more rounded behavior is found when wetting at the boundaries 
occurs. The difference between incompletely wet states and completely wet states is even more pronounced when one 
studies the magnetization fluctuation: for wet states, this fluctuation is very large for almost all values of D until the 
transition in the bulk is reached (Fig. 116b ). This behavior is corroborated by an examination of snapshot pictures 

(Fig.nzD. 

Finally, Fig. [15] shows an attempt to study the divergence of the interfacial adsorption, that we expect when we 
approach the bulk transition in the wet phase in the region where a flrst-order transition occurs. However, it is seen 
that the effective slope of these data changes slightly when the transition is approached, and also depends on \Hi\/ J. 
As a consequence, we must conclude that the asymptotic region where a universal exponent could be estimated has 
not been reached. In view of the strong and slow fluctuations of the interface it is not straightforward to cope with 
this problem. 

IV. CONCLUSIONS 

In this paper, we have presented a Monte Carlo study of the wetting behavior of the two-dimensional Blume-Capel 
model, mapping out the surface critical held Hiw{T, D) where critical wetting occurs for a broad range of parameters 
T and I?, both in the region where the transition of the Blume-Capel model in the bulk is of second order (cf. Fig. 
Hg) and where is of first order (cd. Fig. Hg. We suggest that Hi^iT.D) behaves as Hic{T,D) cx (TcbiD) - Tfl'^ 
for all D < Dt and Hic{T — > Tcb{D)) ~ Hl^{D) for D > Dt, and propose a qualitative phase diagram in the space 
of the three variables surface field (-ffi), temperature (T) and "crystal field" (D), as sketched in figure HM Thus, 
we present evidence for the universality of the surface critical exponent Ai along the critical line of the Blume- 
Capel model. However, the precise crossover behavior in the inmediate vecinity of the tricritical point {Dt,Tcb{Dt)) 
could not yet been studied: first of all, the location of this point needs to be more accurately determined, and also 
more powerful numerical methods (rather than the straightforward Monte Carlo simulation method just using the 
Metropolis algorithm) would be required. 

Already for the task that we did achieve we encountered the necessity to reconsider the finite-size scaling approach 
to the study of wetting transitions. Sec. II. We propose that by using a procedure where the linear dimensions L, M 
of the L X M strip are varied such that the generalized aspect ratio c = L'^/M is kept constant, cf. Eq. it is 

possible to analyze the data in full analogy to the study of a bulk phase transition, where one uses antisymmetric 
surface fields (Hi = —Hl). Then, simply the total magnetization m in the system acts like as order parameter, but 
the appropiate critical exponent is /3 = 0. Using the exactly known results for the Ising model in the square lattice, 
which is the limit of the Blume-Capel model for D — ^ oo, as a test case, our new finite-size scaling approach is nicely 
verified. Note that Hic{T,D oo) and the critical exponents ~ 2, i'± — 1, As = 3, ag — 0, /?s — —1, and 7s = 4 
for critical wetting are known. Then, we propose the corresponding exponents when one treats the transition not as 
a surface free energy singularity of the semi-infinite system in the limit where the bulk field H tends to zero, but 
as a bulk singularity of the L x M system for H ^ in the limit L oo, L"^ /M ~ c ~ constant. These expnents 
are a = — 1, /3 = and 7 = 3. We have shown that these exponents satisfy all expected scaling relations and we 
have verified them from our simulations. Use of this formulation of finite-size scaling for critical wetting in c? = 2 
dimensions is a convenient and useful tool; it would be interesting to study the extension of this method to d = 3 
dimensions where = (logarithmic growth) and this extreme shape anisotropy that our method then requires 
makes the task obviously very much harder. 
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Since the Blume-Capel model differs from the Ising model by the presence of vacancies as a third component, it 
is interesting to ask to what extent vacancies get enriched at the interfaces in this situation of interfaces confined 
by walls. We found that the situation is fully analogous to interfaces between coexisting phases in the bulk: in 
the second-order region the interfacial adsorption is finite, while in the first-order region the predicted divergence of 
interfacial adsorption is found. 

Of course, there are many cases where by changing a parameter the order of a phase transition in the bulk changes 
from second to first order, and one can ask how wetting phenomena are affected. We hope that the present work will 
stimulate further theroretical and experimental studies along such lines. 
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FIG. 1: Schematic description of the system geometry and its state slightly below the wetting transition temperature Tw{Hi) 
such that both and ^± are much larger than the lattice spacing. Coarse-graining the local magnetization on a length scale 
intermediate between the lattice spacing and one is left with one coarse-grained contour (the interface) separating the 
domain with positive magnetization, which was assumed to be the majority domain in the figure, without loss of generality, 
so the interface is still bound to the lower wall, from the domain with negative magnetization. The sign of the magnetization 
of the domains is indicated by double arrows. Note that ^|| oc ^i, and the mean distance of the interface from the nearest 



boundary also is of the same order as ^j.. 
x-direction (p.b.c.) is indicated. 



The choice of linear dimensions L, M, and of the periodic boundary conditions in 
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FIG. 2; (Color online). Monte Carlo estimates for the locations of the phase boundary of the two- dimensional Blume-Capel 
model, in the plane of variables fcsT/J (abscissa) and D/J (ordinate). Full squares denote first-order transitions, the full 
diamond denotes the tricritical point {kBTt/J « 0.609, Dt/J ~ 1.965), and full circles denote second-order transitions. Broken 
straight lines connecting these points are guides to the eye only. The ordered ferromagnetic phase (ord) occurs below the line, 
and the disordered phase (dis) above it, as indicated. Note that the standard two-dimensional Ising model results in the limit 
where D/J tends to minus infinity; so its transition is beyond the scale of the diagram. 
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FIG. 3: (Color online). Plots of the average absolute value of the magnetization (|m|) (a), the magnetization square (b), and 
the cumulant (c) versus temperature, for the choice D/J = — oo, H\/J = 0.70, thus all data shown refer to the standard 
two-dimensional Ising model. The vertical lines indicate the exactly known [l^ location of the wetting transition temperature. 
Curves connecting points are drawn as guides to the eye (also in Figs. |4l [5]). 
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FIG. 4: (Color online). Same as Fig. |3] but for D/ J = 0.0, Hi/ J — 0.55. From the intersection points one can conclude 
kBT-uj{Hi) / J — 1.393 ± 0.004, where the error bar merely reflects the scattering of the intersection points of the measured 
observables. 
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FIG. 5: (Color online). Same as Fig. [3l but for D/ J = 1.75, H\/ J — 0.85. From the intersection points one can conclude 
kBTw(H\) / J — 0.538 ± 0.004, where the error bar merely reflects the scattering of the intersection points of the measured 
observables. 
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FIG. 6; (Color online). Snapshot pictures of the spin configuration of the Blume-Capel model at D/ J = 1.5 and H\/ J — 0.7 
for (L, M) = (18, 288) at three temperatures: ksT/J = 0.445 (left), 0.481 (middle) and 0.518 (right). The sites where S^ = -1 
are shown in black, the sites where Si = are shown in grey (red color online), while sites with Si — +1 are left blank. Note 
that the case ksT/ J = 0.481 is close to the wetting transition. 
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FIG. 7: (Color online), a) Extrapolation of the peak position ftsTmax/J versus 1/L (see Eq. (|44[) ) for the case D/J = 
—oo, Hi/ J = 0.7. The full square indicates the exactly known 193 location of the wetting transition temperature, and the 
dashed line has been drawn in order to guide the eye. Note that /M — 9/8 (Eq. (|45])) is chosen throughout, b) Linear-linear 
plot of fcsTxmax vcrsus L'^, for the same choice as used in a), to demonstrate the power law kBTxmsL^ oc L'^ or (M^^^), see 
Eq.[37l 
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FIG. 8: (Color online). Scaling plot of (|m|)T (a), (m^) (b), and U (c) versus ty/M for the same case D/ J = — oo, Hi/ J = 0.7, 
and various choices of L (for /M = 9/8) as indicated in the figure. The vertical line is a reminder that in this case the 
transition point f = is known exactly. 
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FIG. 9: (Color on line). Same as Fig. [HI but for the case D/J — 1.75, Hi/J — 0.85. The transition point t = was determined 
in Fig. [5] 
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FIG. 10: (Color online), (a) Plots of PL,M{'m) versus m as obtained for various choices of L and M, as indicated, but keeping 
L/^ /M — 9/8 (Eq. (|45l) )and for the case D = — oo. Data obtained for Hi/ J = 0.70 and the exactly known [l^ location of the 
wetting transition temperature, namely kBTw{Hi)/ J ~ 1.6111. (b) As in (a) but for the case D/J — 1.50, and Hi/ J = 0.70. 
The best data collapse of the data is obtained for kBTw[H\) / J = 0.792, in agreement with the intersection points of the 
quantities (|m|)T, {'m?)T and U{T), which are not shown here for the sake of space. 
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FIG. 11: (Color online), (a) Plots of PL,M('m) versus m as obtained for L = 30 and various choices of M, in order to scan 
a wide range of the generalized aspect ratio c = I? jM, as indcated. Data obtained for Dj J = —oo, Hi/J = 0.70 and the 
exactly known [l^ location of the wetting transition temperature, namely kBT-w{H\) / J ~ 1.6111. (b) Intersection points of the 
quantities {m^), and U{T), versus the sample generalised aspect ratio c = /M as obtained for the same parameters as 

in (a). The dashed lines showing a linear extrapolation to c = have been drawn for the sake of comparison. 
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FIG. 12: (Color online), a) Plot of Hic{T) / J , the inverse function of kBT^iHi)/ J versus temperature for a range of values of 
D/J, as indicated in the figure, b) Plot of H\c{T)/ J versus {kBTcb{D)/J - ksT/jy^^ for a range of value D/J, to show that 
the exponent Ai — 1/2, irrespective of D/J, in the region of the second-order transition. In both cases a) and b), the full line 
corresponds to the exact solution for the case D/J = — oo given by Eq. ((2]) [l9| . 
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FIG. 13; (Color online). Plot ol Hic{T)/ J versus D/J for the tricritical temperature [ksT / J = 0.609), and three temperatures 
in the region of ksT j J where the transition in the bulk clearly is first order, ksT j J = 0.3045, ksT j J = 0.406, and ksT j J = 
0.500, respectively, as indicated. Also, two sets of data corresponding to the range of ksT j J where the transition in the bulk 
clearly is second order, ksT j J = 0.8000, and ksT j J — 1.00, respectively, are shown as indicated. Note that dashed-straight 
lines connecting the points are drawn as guide to the eye only. 
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FIG. 14: Schematic description of the wetting behavior of the two-dimensional Blume-Capel model in the space of variables 
Hi, D and T. Complete wetting occurs above the surface of critical wetting transition, which ends at Hi = when the bulk 
transition is continuous, and at the line Hl^{T) in the region beyond the tricritical point {Dt,Tt), where the transition in the 
bulk is of 1°* order. 
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FIG. 15: (Color online). Vacancy concentration profiles (a,c,d) and the magnetization profiles (b) plotted as a function of the 
distance from the wall on which the positive surface field acts, in the case where the majority of the thin strip is in a state of 
negative magnetization in case (a), (b). Data refer to the lattice with linear dimensions (L, M) = (24,288) and the reduced 
temperature fcsT/J = 0.406, surface field \H-i\/ J = 0.91, and several choices of the parameter D/ J , as indicated. 
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FIG. 16: (Color online). Total concentration of vacancies in the system (a) and its fluctuation (b) plotted vs. D/J for 
ksT/J = 0.406 and several choices of the surface field h = \Hi\/J, as indicated. Part(b) includes also the fluctuation of the 
magnetization m for comparison. The system linear dimensions were again chosen as {L, M) — (24, 288). 




FIG. 17: (Color online). Typical snapshot pictures of the spin configurations obtained for ksT/ J = 0.406, \Hi\/ J = 0.91 and 
three values of D/J : D/J — 1.80 (left panel, refers to an uncompletely wet state), 1.90 (medium panel, wet phase), 2.02 (right 
panel, disordered phase beyond the first order transition of the bulk). Sites i with 5i = 1 are shown by black dots, sites i with 
Si = are shown as grey (b/w) or red (in color) dots, sites i with Si — —1 are left blank. 
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FIG. 18: (Color online). Log-log plot of the excess density of the vacancies versus (-Dcoex — D)/J at ksT/J = 0.406 (where 
Dcoc^/J = 1.996) and different values of the surface field h = H\/J. Broken lines show the theoretical slopes 1/3 (expected 
in the first-order region far off from the tricritical point) and 4/9 (expected at the tricritical point), see also Eq. ((Oj. Chosen 
lattice size was {L,M) = (201,400). 



29 



[1] J.S. Rowlinson and B. Widom, Molecular Theory of Capillarity (Oxford University Press, Oxford, 1982) 
[2] P.G. De Gennes, F. Brochard-Wyart, and D. Quere: Capillarity and Wetting Phenomena: Drops, Bubbles, Pearls, Waves 
(Springer, Berlin, 2003) 

[3] M. Schon and S. Klapp, Nanoconfined Fluids: Soft Matter Between Two and Three Dimensions (J. Wiley & Sons, New 
York, 2006) 

[4] I. Brovchenko and A. Oleinikova, Interfacial and Confined Water (Elsevier, Amsterdam, 2008) 

[5] G. Docker and J.B. Schlenoff (eds.) Multilayer Thin Films: Sequential Assembly of Nanocomposite Materials (Wiley- VCH, 
Wcinheim, 2002) 

[6] Y. Champion and H.-J. Fccht (eds.) Nano-Architectured and Nano-Structured Materials (Wiley- VCH, Wcinheim, 2004) 

[7] S. Zhang (ed.) Handbook of Nano structured Thin Films and Coatings, Vol. 1-3 (CRC Press, Boca Rotan, 2010) 

[8] T.M. Squires and S.R. Quake, Rev. Mod. Phys. 77, 977 (2005) 

[9] P.G. De Gcnncs, Rev. Mod. Phys. 57, 827 (1985) 
[10] D.E. Sulhvan and M.M. Telo da Gama, in C. Croxton (ed) Fluid Interfacial Phenomena (Wiley, New York, 1986) p. 45 
[11] D.B. Abraham, in C. Domb and J.L. Lebowitz (eds.) Phase Transitions and Critical Phenomena, Vol. 10 (Academic, 
London, 1986) p.l 

[12] S. Dietrich, in C. Domb and J.L. Lebowitz (eds.) Phase Transitions and Critical Phenomena, Vol. 12 (Academic, London, 
1988) p.l. 

[13] M. Schick, in J. Charvolin, J.-F. Joanny, and J. Zinn-Justin (eds.) Liquids at Interfaces (Elsevier, Amsterdam, 1990) p. 
415 

[14] G. Forgacs, R. Lipowsky and T.M. Niouwcnhuizon, in C. Domb and J.L. Lebowitz (eds.) Phase Transitions and Critical 

Phenornana, Vol. 14 (Academic, London, 1991) Chap. 2 
[15] D. Bonn and D. Ross, Rep. Progr. Phys. 64, 1085 (2001) 

[16] D.R. Clarke, M. Riihle, A. P. Tomsia (eds.) Low and High Temperature Wetting: State of the Art; Ann. Revs. Mater. Res., 

Vol. 38 (Ann. Revs., Palo Alto, 2008) 
[17] D. Bonn, J. Eggcrs, J. Indckcu, J. Mcunicr, and E. RoUey, Rev. Mod. Phys. 81, 739 (2009) 
[18] K. Binder, D.P. Landau, and M. Miillcr, J. Stat. Phys. 110, 1411 (2003) 
[19] D.B. Abraham, Phys. Rev. Lett. 44, 1165 (1980) 
[20] M.E. Fisher and H. Nakanishi, J. Chem. Phys. 75, 5857 (1981) 
[21] H. Nakanishi and M.E. Fisher, Phys. Rev. Lett. 49, 1565 (1982) 
[22] H. Nakanishi and M.E. Fisher, J. Chem. Phys. 78, 3270 (1983) 
[23] K. Binder, D.P. Landau and D.M. KroU, Phys. Rev. Lett. 56 , 2272 (1986) 
[24] D.B. Abraham and E.R. Smith, J. Stat. Phys. 43, 621 (1986) 
[25] K. Binder and D.P. Landau, Phys. Rev. B37, 1745 (1988) 
[26] K. Binder, D.P. Landau, and S. Wansleben, Phys. Rev. B40, 6971 (1989) 
[27] E.V. Albano, K. Binder, D.W. Heermann, and W. Paul, Surf. Sci. 223, 151 (1989) 
[28] E.V. Albano, K. Binder, D.W. Heermann, and W. Paul, J. Stat. Phys. 61, 161 (1990) 
[29] A.O. Parry and R. Evans, Phys. Rev. Lett. 64, 439(1990) 
[30] A.O. Parry and R. Evans, Physica A 181, 250 (1992) 

[31] K. Binder, D.P. Landau, and A.M. Forrenberg, Phys. Rev. Lett. 74, 298 (1995) 

[32] K. Binder, D.P. Landau, and A.M. Forrenberg, Phys. Rev. E51, 2823 (1995) 

[33] K. Binder, R. Evans, D.P. Landau, and A.M. Ferrenberg, Phys. Rev. E53, 5023 (1996) 

[34] A. Maciolek and J. Stecki, Phys. Rev. B54, 1128 (1996) 

[35] A. Maciolek, J. Phys. A29, 3837 (1996) 

[36] E. Carlon and A. Drzewinski, Phys. Rev. E57, 2626 (1998) 

[37] E.V. Albano, K. Binder, and W. Paul, J. Phys.:Condens. Matter 12, 2701 (2000) 
[38] A. De VirgiUis, E.V. Albano, M. Miiller, and K. Binder, Physics A352, 477 (2005) 
[39] B.J. Schulz, K. Binder, and M. Miiller, Phys. Rev. E71, 046705 (2005) 
[40] L. Pang, D.P. Landau, and K. Binder, Phys. Rev. Lett. 106, 236102 (2011) 

[41] L.D. Gelb, K.E. Gubbins, R. Radhakrishnan, and M. Slivinska-Bartkowiak, Rep. Progr. Phys. 63, 1573 (1999) 
[42] E.V. Albano, K. Binder, D.W. Heermann, and W. Paul, J. Chem. Phys. 91, 3700 (1989) 
[43] K. Binder and D.P. Landau, J. Chem. Phys. 96, 1444 (1992) 

[44] O. Dillmann, W. Janke, M. Miiller, and K. Binder, J. Chem. Phys. 114, 5823 (2001) 
[45] K. Binder and P.O. Hohenberg, Phys. Rev. B6, 3461 (1972) 
[46] K. Binder and P.C. Hohenberg, Phys. Rev. B9, 2194 (1974) 

[47] K. Binder, in C. Domb and J.L. Lebowitz (eds.) Phase Transitions and Critical Phenomena, Vol. 8 (Academic, London, 
1983) p.l. 

[48] L. Onsager, Phys. Rev. 65, 117 (1944) 
[49] M. Blume, Phys. Rev. 141, 517 (1966) 
[50] H.W. Capel, Physica 32, 966 (1966) 

[51] S. Sarbach and I.D. Lawrie, in C. Domb and J.L. Lebowitz (eds.) Phase Transitions and Critical Phenomena, Vol. 9 
(Academic, London, 1984) p.l 



30 



[52] K. Binder and D.P. Landau, Surf. Sci. 61, 577 (1976) 

[53] J.L. Cardy, in C. Dorrib and J.L. Lebowitz (eds.) Phase Transitions and Critical Phenomena, Vo. 11 (Academic, London, 
1987) p.55 

[54] R.B. Pearson, Phys. Rev. B22, 2579 (1980) 

[55] B. Nienhuis, J. Phys. A: Math. Gen. 15, 199 (1982) 

[56] H.E. Stanley, An Introduction to Phase Transitions and Critical Phenomena (Oxford Univ. Press, Oxford, 1971) 

[57] K. Binder, in D.G. Pettifor (ed.) Cohesion and Structure of Surfaces (Elsevier, Amsterdam, 1995) p. 121. 

[58] W. Seiko and J. Ycomans, J. Phys. A16, 2789 (1983) 

[59] W. Sell^e, D.A. Huso, and D.M. Kroll, J. Phys. A17, 3019 (1984) 

[60] W. Sellie, Surf. Sci. 144, 176 (1984) 

[61] R. Lipowsky, Phys. Rev. B32, 1731 (1985) 

[62] D.M. Kroll, R. Lipowsky, and R.K.P. Zia, Phys. Rev. B32, 1862 (1985) 

[63] A.O. Parry, J.M. Romero-Enrique and A. Lazarides, Phys. Rev. Lett. 93, 086104 (2004) 

[64] A.O. Parry, C. Rascon, N.R. Bernardino, and J.M. Romero-Enrique, J. Phys.: Condens. Matter 18, 6433 (2006); ihid 19, 
416105 (2007) 

[65] A.O. Parry, C. Rascon, N.R. Bernardino, and J.M. Romero-Enrique, Phys. Rev. Lett. 100, 136105 (2008) 

[66] A.O. Parry, J.M. Romero-Enrique, N.R. Bernardiono, and J. Rascon, J. Phys.: Condens. Matter 20, 505102 (2008) 

[67] K. Binder and J.S. Wang, J. Stat. Phys. 55, 87 (1989) 

[68] K. Binder, in V. Privman (cd.) Finite Size Scaling and Numerical Simulation of Statistical Systems (World Scientific, 

Singapore, 1990) p. 173 
[69] S.M. Bhattacharjee and J.F. Nagle. Phys. Rev. A, 31, 3199 (1985) 
[70] S.M. Bhattacharjee and J.J. Rajasekaran. Phys. Rev. A, 44, 6202 (1991) 
[71] K. Binder, Phys. Rev. Lett. 47, 693 (1981) 
[72] K. Binder, Rep. Progr. Phys. 60, 487 (1997) 

[73] L.D.C. Jaubert , J.T. Chalker, P.C.W. Holdsworth, and R. Moessner. Phys. Rev. Lett. 105, 087201 (2010) 
[74] A. Milchev, M. Miiller, K. Binder, and D.P. Landau, Phys. Rev. Lett. 90, 136101 (2003) 
[75] A. Milchev, M. Miiller, K. Binder, and D.P. Landau, Phys. Rev. E68, 031601 (2003) 
[76] A. Milchev, M. Miiller, K. Binder, Europhys. Lett. 70, 348 (2005) 
[77] A. Milchev, M. Miiller, and K. Binder, Phys. Rev. E72, 031603 (2005) 

[78] A. DeVirgihis, E.V. Albano, M. Miiller, and K. Binder, J. Phys.: Condens. Matter 17, 4579 (2005) 

[79] A. Winkler, D. Wilms, P. Virnau, and K. Binder, J. Chem. Phys. 133, 164702 (2010) 

[80] D.P. Landau and R.H. Swendscn, Phys. Rev. B33, 7700 (1986) 

[81] J.D. Kimel, S. Black, P. Carter, and Y.L. Wang, Phys. Rev. B35, 3347 (1987) 

[82] W. Speth, Z. Phys. B51, 361 (1983) 



